Nematic domains and resistivity in an itinerant metamagnet coupled to a lattice 
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The nature of the emergent phase near a putative quantum critical point in the bilayer ruthenate 
Sr3Ru207 has been a recent subject of intensive research. It has been suggested that this phase 
may possess electronic nematic order(ENO). In this work, we investigate the possibility of nematic 
domain formation in the emergent phase, using a phenomenological model of electrons with ENO 
and its coupling to lattice degrees of freedom. The resistivity due to the scattering off the domain 
walls is shown to closely follow the ENO parameter. Our results provide qualitative explanations 
for the dependence of the resistivity on external magnetic fields in Sr3Ru207. 

PACS numbers: 71.27.+a,75.30.Kz,71.10.-w,63.20.Kr 



Quantum critical behavior in metamagnetism repre- 
sents the most unusual critical phenomenon in itinerant 
electron systems. It was shown in a series of remark- 
able experiments that the critical end point associated 
with the metamagnetic transition(MMT) can be driven 
to zero temperature in the bilayer ruthenate Sr3Ru207 by 
changing the angle between the external magnetic field 
and the c-axis of the material. [1, 2] The transport and 
thermodynamic measurements reveal unusual tempera- 
ture dependence near the putative quantum critical end 
point, indicating that the underlying itinerant electrons 
are strongly interacting with the quantum critical fluctu- 
ations of the magnetization. [1-4] 

Later experiments on the ultra-pure samples, however, 
have shown that the quantum critical end point is eventu- 
ally avoided when the finite temperature critical point is 
pushed further down to zero temperature, and is replaced 
by two consecutive MMTs.[5-7] The phase bounded by 
two MMTs exhibits unusually high resistivity and its na- 
ture has been a recent subject of intensive research. 

Currently there exist two competing proposals for the 
nature of this emergent phase. Binz et al. [8] proposed a 
phenomenological Ginzburg-Landau theory of magnetic 
Condon domains and suggested that the high resistivity 
may come from the scattering of electrons off the mag- 
netic domains. Another proposal involves the formation 
of ENO in the emergent phase; namely the Fermi sur- 
face of the electrons spontaneously breaks the lattice ro- 
tational symmetry. [5, 9] Two possible orientations of the 
ENO may lead to domain formation and may indeed lead 
to high resistivity due to domain wall scattering. 

Very recently it has been reported that the high re- 
sistivity seen in the emergent phase is insensitive to de- 
magnetization factor or shapes of the samples. [10] This 
is at odds with a key prediction of the magnetic Condon 
domain theory; the system is supposed to be very sen- 
sitive to demagnetization factor. [8] This may leave only 
the theory of ENO as an alternative proposal. 

In this paper, we offer a theory of the nematic domain 
formation and the resulting high resistivity in the emer- 



gent phase. The possibility of the nematic domains has 
not been studied theoretically even though such possi- 
bility was briefly discussed in the literature. [5, 9] Our 
objective in the current work is to provide a minimal 
theoretical description that can explain key experimen- 
tal observations at qualitatively level. 

In search for a concise description, we start with an 
effective electronic Hamiltonian that captures the essen- 
tial physics of ENO, namely we consider the quadrupolar 
interaction between electrons. [9, 11] This interaction can 
be regarded as the angular-momentum-two channel of 
some general interaction or the Landau quasiparticle in- 
teraction. In real materials, the interaction in all angular 
momentum channels may exist, not just in the quadrupo- 
lar channel. The dominance of the quadrupolar channel 
is an assumption of our phenomenological theory and its 
validity should be tested by comparing the results with 
the experimental data. While it has not been proven, 
it is conceivable that such an instability or a dominant 
interaction is more likely to happen near an itinerant 
quantum critical point because the Fermi surface would 
be very soft due to strong critical fluctuations. 

In view of the magnetostriction data indicating a 
close relationship between the lattice and the emergent 
phase [5], we also consider the electron-lattice interac- 
tion. Since we are studying an effective Hamiltonian, we 
are only interested in general properties of this Hamil- 
tonian and the precise values of the parameters should 
not be taken seriously. Remarkably this simplified model 
contains essential ingredients of the nematic domain for- 
mation and the resulting high resistivity. It is found that 
the electron-lattice interaction makes the formation of 
the nematic domains (at finite temperatures) much more 
likely by offering a large number of metastable domain 
configurations. Another key result of our work is that 
the magnetic field dependence of the resistivity in the 
nematic domain phase closely follows the field depen- 
dence of the ENO parameter. Thus the resistivity can 
be regarded as a measurement of the ENO. 

Mean field theory of the nematic order: The 



simplest model for the ENO on the square lattice can be 
written by keeping only the interactions in the quadrupo- 
lar channel [11] or the lattice equivalent of the angular 
momentum I = 2 channel in the continuuni[12]. 
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where ^k = —t(coskx + cos ky) — 2t' cosk^ cos ky — /i is 
the single-particle dispersion and a = 1, 2. /ki = cosfc^: — 
cos ky and /k2 = sin k^ sin ky are the form factors in the 
quadrupolar channels. F{q) represents a short range in- 
teraction in real space with F(q ^0) = F = constant. 
Here we suppress the spin indices and it should be un- 
derstood that the interaction terms exist in both spin 
channels with possibly different interaction strength. 

The ENO parameter is Q(q) = -F^k /ki(Ck+qCk)- In 
the unform state q -^ 0, this order parameter leads to the 
renormalized dispersion in the mean field theory, ^^""^ = 
— {t + Q{0)) cos kx — (i — (5(0)) cos ky — 2t' cos kx cos ky — fi. 
Clearly, the resulting Fermi surface would break the 7r/2 
rotation symmetry of the lattice. An alternative ENO pa- 
rameter with /k2 vanishes for the interaction form chosen 
here. [12] 

In the presence of a magnetic field, the Zeeman energy 
leads to different chemical potentials for the spin-up and 
-down electrons, /i|^| = M ± Ms^, where hb is the Bohr 
magneton and H is the external magnetic field. As the 
magnetic field increases, the Fermi surface volume of the 
majority (minority) spin increases (decreases). When the 
Fermi surface of the majority or minority spin touches 
the van-Hove singularity, a first order phase transition 
to a nematic phase occurs with a sudden opening and 
elongation of Fermi surface along a- or 6-axis.[9] Here 
the Fermi surface volume of the majority or minority 
spin jumps, leading to an MMT. When the magnetic field 
increases further, the other side of Fermi surface touches 
the van-Hove singularity again. This leads to another 
first order transition and the lattice rotational symmetry 
is restored afterwards. [9] 

In the uniform nematic state, the spontaneously cho- 
sen nematic direction can be either the a- or b- axis. On 
the other hand, non-uniform solutions can arise as meta- 
stable states in the mean field theory. In particular, the 
nematic domains with different orientations can coexist. 
These static domain structures, however, may not be fa- 
vorable even as meta-stable states in two and one dimen- 
sions because of large quantum fluctuations. Thus if the 
nematic domains were to arise, there should be additional 
degrees of freedom. 

Normal modes of lattice distortions: Magne- 
tostriction effect observed in Sr3Ru2 07 suggests that the 
electron-lattice interaction is a natural degree of freedom 
to consider. In the uniform nematic-ordered state, the 
electron-lattice interaction would induce a tetragonal lat- 
tice distortion. This system-wide shape change, however. 



is energetically very costly and the system may rather 
prefer to have twin boundaries between the phases with 
different nematic orientations so that the overall shape 
is preserved. Thus the electron-lattice interaction may 
indeed prefer the formation of the nematic domains. 
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FIG. 1: Normal modes 

We treat the lattice classically assuming slow lattice 
dynamics. For simplicity, we will only consider two- 
dimensional lattice distortions even though the full anal- 
ysis of the electron-lattice coupling may require a three 
dimensional lattice structure. Implications to the distor- 
tions out of plane will be discussed later. 

It is useful to write the lattice Hamiltonian in terms of 
distortion modes of a lattice plaquette. Apart from the 
rotational and two translational modes of the entire pla- 
quette, there are five normal modes described in Fig. 1. 
These normal modes are, however, not independent in the 
lattice and three constraints between the normal modes 
are necessary. [13] Using the classical approximation, we 
can write the harmonic part of the lattice Hamiltonian 
as follows. 
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where it should be understood that only two of the nor- 
mal modes are independent. 

It is easy to have a grip on the physical meaning of the 
parameters in this Hamiltonian. The normal modes s± 
appear at the domain walls so that B is related to the 
domain wall energy. It is also clear that 63 mode cou- 
ples directly to the ENO; the domain boundaries of the 
lattice distortions related to the 63 modes will coincide 
with those of the ENO. It is expected that the 63 mode 
is particularly affected and will be soft. Thus we assume 
that A3 is small while all the An and B are positive. 

Using the constraints, the lattice Hamiltonian with an 
anharmonic correction can be rewritten as follows. 



^lat — 



iiat - E 

q, a, 6— 1,3 



ea(q) — ;; — 66 (q)- 



c. 



Y.^eM\ (3) 



/3a /3b 



B 



Pi 



p2[l3l6ab + I33{1 - Sab)]- 

1 - COS qx COS qy, /32(q) = 
cosgj^, and Pi{q) = (1 — 



with Mab = Aa5ab + A2 

C3 is a constant and /3i(q) 
- sin qx sin qy, /33(q) = cos 9. 
cosqx){l -cosqy). 

Electron-lattice interaction: We consider the Su- 
Schrieffer-Heeger type electron-lattice interaction, 

i^ci-iat = 5 E[(^^+* " ^")cl+£C» + {x^y)+ h.c], (4) 



where di represents the lattice displacement. Writing this 
interaction in terms of the lattice normal modes, we get 
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where /ko = cos k^ + cos ky . Notice that the 63 mode 
couples to the ENO through the form factor /ki and the 
ei mode to the kinetic energy. 

Domain structures: We numerically determine the 
electronic and lattice configurations of the system de- 
scribed hy H = Hc\ + -ffci-iat + ^lat • The simplest domain 
wall solutions can be obtained by assuming that the spa- 
tial modulations exist only along the diagonal direction. 
This makes the numerical works easier while one can still 
understand the key properties of the domain wall config- 
urations. We also assume that the ei mode is much stiffer 
than the 63 mode so that ei mode is not heavily affected 
by the nematic transitions. We will discuss later the role 
of the ei mode in the context of the magnetostriction 
measurements in Sr3Ru207. 

The non-uniform ENO parameter and lattice configu- 
rations are determined by the self-consistent mean field 
theory described below. We compute the energy of the 
anharmonic part of the lattice Hamiltonian in real space 
while the rest of the energy is computed in the reciprocal 
space. We use the periodic boundary condition in the 
reciprocal space and the fast-Fourier transform to con- 
nect the real and reciprocal space computations. Here 
the number of discrete q points in the Brillouin zone of 
the reciprocal space is related to the lattice size in the 
real space. We first consider Hd + -ffei-iat and determine 
the ENO parameter Q(q;e3) of the electronic state for 
a given initial lattice configuration {63}. The resulting 
electronic state |<i>(e3))ci is used to compute the energy 
^01(63) = ci($(e3)|-ffoi + -ffci-iat(e3)|$(e3))ci. Next we 
minimize the total energy i?totai(e3) = Ec\{e-i) + H\a.t{e?,) 
by the Eulcr method[13]. The resulting lattice config- 
uration is in general different from the initial one and 
it is then fed back to determine a new electronic state 
$(63)). We update the lattice configuration {63} fol- 
lowing the maximum energy slopes with respect to 63 (q) 
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FIG. 2: (Color online) The red, blue, and green represent 
the negative, positive, and zero ENO parameter respectively. 
Here, Ai = 20, A2 = 20, A3 = 0.001, B = 1, C3 = 100, 
t = 0.1, t' = Q, F = 0.2, and g = 0.2. The electron filling is 
il§ per site to be slightly away from the van Hove singularity. 



The relative probability of getting the metastable do- 
main solutions is obtained by trying a large number (typ- 
ically 50 in a given set) of random initial conditions. This 
probability is only 32% for the 8x8 lattice, it increases 
to 50% for the 18 x 18 lattice, and 63% for the 20 x 20 lat- 
tice with 2^3 % error. Thus the number of metastable 
domain solutions increases as the system size increases 
in accord with similar computations in manganites.[13] 
These domain solutions imply a large number of local 
minima in the free energy and will be likely chosen at 
finite temperatures by entropic effect. Fig. 2 shows a 
typical uniform and a meta-stable domain solution. The 
qualitative features of these solutions do not depend on 
the details of the parameters chosen here; for example, 
the electron filling factor Iff with i' = is chosen just 
to avoid the van Hove singularity and a similar result is 
obtained for the half-filling with a small finite t' . 

The stability of domains: The stability of the do- 
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FIG. 3: (Color online) k as a function of (a) the electron- 
lattice coupling and (b) the anharmonicity. We use t' — 0, 
F = 2t, and B = 20t. C3 = 0.05 for (a) and g = 0.1 for (b). 

main configurations would be determined by the barrier 
potential in the energy landscape and the energy differ- 
ence between the uniform state and the domain structure. 

Here we consider an indirect criterion of the 'local' sta- 
bility, K = ^ -g^A^Etotai evaluated at the local min- 
ima corresponding to the meta-stable domain configura- 
tions. This quantity is the curvature of the energy topog- 
raphy in the hyperspace of lattice configurations. Larger 
K would correspond to a deeper local minimum, implying 
perhaps a more stable configuration. Even though this 
quantity measures the stability of the domains in the 
lattice deformations, the nematic domains would follow 
because they are basically locked together. 

Of particular interest is the stability with respect to the 
electron-lattice coupling g and the strength of the anhar- 
monicity C3. Fig. 3 (a) and (b) show the dependence of 
K on the dimensionless coupling constants g = ge^/t and 
C3 — C^e^/t, where 63 is the average of 63 and about the 
order of 0.1 lattice constant. The increase of k in both 
of the uniform and domain solutions may lead to bigger 



barrier height and the faster increase of k for the domain 
configuration may imply that the domain structure at 
large g and C^ would be at a relatively deeper minimum. 
Resistivity: The large enhancement of the residual 
resistivity in the emergent phase may be explained by 
the scattering at the nematic domain walls. We adopt the 
Landauer formula to estimate the conductance through 

2 

a domain wall, namely Gdw = ^T, where T is the 
transmission coefficient. [14] In the case of a domain wall 
parallel to the 6-axis, the transmission coefficient along 
the a-axis (or the x-axis) is given by T = X^ks-FS-^k, 
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and k e FS ensures that only the k vectors residing 
on the appropriate Fermi surface are included. Here 
^p '{k) = -§^i^ 'ikGFS, Vq is the strength of the scat- 
tering potential at the interface, and ± sign denotes the 
direction of the nematic-ordered Fermi surfaces in the dif- 
ferent sides of the domain wall. If the scattering processes 
at different domain walls are incoherent, the resistivity 
due to A^DW domain walls would be 

Pdw oc iVDwG;^^ = A^DW /(uo - uiQ^) (6) 



where wq 
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sistivity is shown in Fig. 4 for i' ^ and t' = cases. 
Notice that i' ^ case leads to an asymmetric profile 
of the resistivity similar to that seen in the experiment. 
The important point, however, is that the overall shape 
of the resistivity follows the behavior of the ENO and 
this aspect is independent of the parameters. 
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FIG. 4: (Color online) The dashed line denotes the uniform 
ENO parameter in the mean-field theory. [12] The solid lines 
represent the resistivity with i' = 0, and 0.4f. Here we use 
p/p{Q = 0) = ^_ " p a with a = 0.24 and F = 1.974f. 

Summery and Discussion: The theory of ENO is 
consistent with several key experimental observations in 
Sr3Ru207.[5, 10] 1) There exist two first order MMTs as 
the magnetic field increases, 2) the electrons are itinerant 
in the entire phase diagram, 3) the emergent phase exists 
only in a window of magnetic fields. 4) the finite tem- 
perature transition to the emergent phase seems to be 



continuous. [5, 9] In this paper, we have shown that this 
theory combined with the coupling to the lattice leads 
to nematic domain formation and this naturally explains 
the high resistivity in the emergent phase. [5, 10] In par- 
ticular, it is shown that the resistivity is determined by 
the magnetic field dependence of the ENO parameter. 

Some additional remarks are in order. The magne- 
tostriction data on Sr3Ru207 show a sudden change in 
the c-axis lattice constant upon entering and exiting the 
emergent phase. [5] This may be explained by a relatively 
soft ei mode that is related to a change in the a6-plane 
unit-cell area. Such a change would lead to elongation of 
the lattice along the c-axis to conserve the three dimen- 
sional unit-cell volume. The jump in the kinetic energy 
of the electrons with the majority/minority spin, via the 
coupling to the ei mode, would induce an abrupt change 
of the c-axis lattice constant at the nematic transitions. 

A very recent experiment observed an interesting de- 
pendence of the resistivity on the relative orientation of 
the in-plane magnetic field and the current direction. [10] 
When they are parallel (perpendicular) to each other, 
the high resistivity disappears (is maintained). This is 
reminiscent of the situation in the nematic quantum Hall 
states in high Landau levels. [15] A similar mechanism [16] 
of the nematic domain alignment by an in-plane mag- 
netic field may work here as well. This similarity, there- 
fore, may make the relevance of the nematic theory even 
stronger. It would be an excellent subject of future study. 
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